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ABSTRACT 

We have investigated the time variations in the light curves from a sample of long and 
short Fermi/GBM Gamma ray bursts (GRBs) using an impartial wavelet analysis. 
The results indicate that in the source frame, the variability time scales for long bursts 
differ from that for short bursts, that variabilities on the order of a few milliseconds 
are not uncommon, and that an intriguing relationship exists between the minimum 
variability time and the burst duration. 
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1 INTRODUCTION 

The prompt emission from Gamma-ray Bursts (GRBs) 
shows complex time profiles that have eluded a gener ally 
accepted explanation. iFenimore fc Ramirez-Ruia (|2000l ) re- 
ported a correlation between variability of GRBs and the 
peak isotropic luminosity. The existence of the variability- 
luminosity correlation suggests that the prompt emission 
light curves have embedded temporal information related 
to the microphysics of GRBs. Several models have been 
proposed to explain the observed temporal variability of 
GRB li ght curves. Leading models such as the internal shock 
model dKobavashi. Piran. and Saril 1 19971 ) and the photo- 
spheric model (|Rvde!l2004 ) link the rapid variability directly 
to the activity of the central engine. Others invoke rela- 
tivistic outflow mechanisms to suggest that local turbulence 
amplified through Lorentz boosting leads to causally dis- 
connected regions which in turn act as independent centers 
for the o bserved prompt emiss ion. Within more recent m od- 
els, both lMorsonv et all (|2010l ) and lZhang fc Yanl |201ll ') ar- 
gue that the temporal variability may show two different 
scales depending on the physical mechanisms generating the 
prompt emission. 

In order to further our understanding of the prompt 
emission phase of GRBs and to explicitly test some of the 
key ingredients in the various models, it is clearly important 
to extract the variability for both short and long gamma-ray 
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bursts in a robust and unbiased manner. It is also clear that 
the chosen methodology should not only be mathematically 
rigorous but also be sufficiently flexible to apply to transient 
phenomena with multiple time scales and a wide dynamic 
range. 

In this paper, we extract variability time scales for 
GRBs using a method based on wavelets. The technique for 
such a temporal analysis is universal, and has the advantage 
over Fourier analysis that transients and frequency correla- 
tions can be more easily picked out in the data. The layout 
of the paper is as follows: the source of the data is described 
in section 2; the main aspects of the wavelet methodology 
are outlined in section 3; in section 4 we provide the details 
of the data analysis; in section 5 we present and discuss our 
main findings; finally, in section 6, we summarize our con- 
clusions. 



2 DATA 

The Gamma-Ray Burst Monitor (GBM) on board Fermi 
observes GRBs in the energy range 8 keV to 40 MeV. The 
GBM is composed of 12 thallium-activated sodium iodide 
(Nal) scintillation detectors (12.7 cm in diameter by 1.27 
cm thick) that are sensitive to energies in the range of 8 keV 
to 1 MeV, and two bismuth germanate (BGO) scintillation 
detectors (12.7 cm diameter by 12.7 cm thick) with energy 
coverage between 200 keV and 40 MeV. The GBM detectors 
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Figure 1. GBM GRB080925775. Preburst portion of the light 
curve, used for background removal, is shown in gray. The burst 
portion, from which a time scale is extracted, is shown in black. 

are arranged in such a way that th ey provide a significant 
view of the sky (|Meegan et al.ll2009f ). 

In this work, we have extracted light curves for the 
GBM Nal detectors over the entire energy range (8 keV 
- 1 MeV, also including the overflow beyond 1 MeV). Typ- 
ically, the brightest three Nal detectors were chosen for the 
extraction. Lightcurves for both long and short GRBs were 
extracted at a time binning of 200 microseconds. The long 
GRBs were extracted over a duration starting from 20 sec- 
onds before the trigger and up to about 50 seconds after 
the Tgo for the burst without any background subtraction. 
For short GRBs, durations were chosen to be 20 seconds 
before the trigger and 10 seconds after the T90. The T90 
dur ations were obtained from the Fermi GBM-Burst Cata- 
log jPaciesas et alJl2012T ). 



3 METHODOLOGY 

We report on our model-independent statistical investiga- 
tion of the variability of Fermi/GBM long and short GRBs. 
We extract this information by using a fast wavelet trans- 
formation to encode GRB light curves into a wavelet rep- 
resentation and then compute a statistical measure of the 
variance of wavelet coefficients over multiple time-scales. 

3.1 Minimum Variability Time Scales 

It is often the case when multiple processes are present that 
one process will dominate the others at certain time scales 
but those same processes may exchange dominance at other 
time scales. A wavelet technique is useful in these situations 
because the variances of wavelet coefficients are sensitive 
to whichever processes dominate the light curve at a given 
time scale. Moreover, the technique can be used to clas- 
sify those dominant processes as well as provide a means 
to determine the characteristic time scale, 773, for which the 
processes exchange dominance. Determination of Tp helps in 
the development of theoretical models and the understand- 
ing of observational data. Indeed, if there is a transition from 
a time-scaling region to that of white noise then there is a 
smallest variability time for the physical processes present. 



3.2 Wavelet Transforms 

Wavelet transformations have been shown to be a natu- 
ral too l for multi-reso l ution analysis of no n-stationary time- 
series (|Flandrinlll989l . Il992l ; lMalIat|ll989h . Wavelet analysis 
is similar to Fourier analysis in many respects but differs in 
that wavelet basis functions are well-localized, i. e. have com- 
pact support, while Fourier basis functions are global. Com- 
pact support means that outside some finite range the am- 
plitude of wavele t basis function s goes to zero or is otherwise 
negligibly small (|Percivalll2000l ). In principle, a wavelet ex- 
pansion forms a faithful representation of the original data, 
in that the basis set is orthonormal and complete. 



3.2.1 Discrete Dyadic Wavelet Transforms 

Given the discrete nature of the data, we employ a dis- 
crete wavelet analysis. The rescaled-translated nature of 
the wavelet basis functions make the wavelet transform 
well-localized in both frequency and time, which results in 
an insensitivity to background photon counts expressed by 
polynomial fits. The level of insensitivity, formally known 
as the vanishing moment condition, can be adjusted by 
the choice of wavelet basis function. By construction, the 
discrete wave let transform is a multi-resolution operation 
(Mallat 1989). Such wavelets, denoted ipj,k{t), form a dyadic 
basis set, i.e. wavelets in the set have variable widths and 
variable central time positions. 

The wavelet analysis employed in this study, as with 
the fast Fourier transform, begins with a light curve with N 
elements, 

Xi = {X ...X N -i}, (1) 

where N is an integer power of two. The light curve is con- 
volved with a scaling function, 4>j,k{U), and wavelet func- 
tion, ipj,k{ti) which are rescaled and translated versions of 
the original scaling and wavelet functions 4>(U) = <f>o,o, and 
ip(U) = V'o.o- Translation is indexed by k and rescaling is 
indexed by j. The rescaling and translation relation is given 
by 

iij, k (t) = 2- j/2 ij(2- j t-k). (2) 

The precise forms of the scaling and wavelet functions 
are not unique. The choices are made acco r ding to the fea - 
tures one wishes to exploit |Percivall |2000| ; lAddisor] 120021 ') . 
The scaling function acts as a smoothing filter for the input 
time-series and the wavelet function probes the time-series 
for detail information at some time scale, At, which is twice 
that of the finest binning of the data, Tbi n . In the analysis, 
the time scale is doubled 

At 2At 

and the transform is repeated until 
At = NT Ua . 

In this analysis we choose the Haar (|Addisonl [2051 ) 
scaling/wavelet basis because it has the smallest possible 
support, has one v anishing moment, and is equivalent to 
the Allan variance (|Howe fc Percival Il995l ), allowing for a 
straightforward interpretation. 
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3.2.2 The Haar Wavelet Basis 

Convolving the light curve, X, with the scaling functions 
yields approximation coefficients, 



a.j,k = {(j)j.k,X). 



(3) 



Interrogating X with the wavelet basis functions yields 
scale and position dependent detail coefficients, 



dj,k = (ipj,k,X), 



(4) 



It is interesting to point out that for the trivial N = 2 case 
the Haar wavelet transform and the Fourier transform are 
identical. 



3.3 Logscale Diagrams and Scaling 

Logscale diagrams are useful for identifying scaling and noise 
regions. Construction of a logscale diagram for e ach GRB 
proce eds from the variance of detail coefficients (|Flandrinl 



(5) 



where the nj are the number of detail coefficients at a par- 
ticular scale, j. A plot of log 2 variances versus scale, j, takes 
the general form 



log 2 Pj = aj + constant, 



(6) 



and is known as a logscale diagram. A linear regression is 
made of each logscale diagram and the slope parameter, 
a, (depicting a measure of scaling) is estimated. White- 
noise processes appear in logscale diagrams as flat regions 
while non-stationary processes appear as sloped regions 
wi th the following condition on the slope param eter, a > 
1 jAbrv et al.ll2003l ; |Perciva3l200d : iFlandrinll 19891 ). 



4 DATA ANALYSIS 

4.1 Background Subtraction 

We now present a method for removing photometric back- 
ground due to noise not intrinsic to the GRB so that phys- 
ical variability arising from the GRB remains for further 
analysis. Background subtraction for a statistical analysis 
of variability via wavelet transforms should proceed in the 
space variances as opposed to a traditional flat or linear 
subtraction of counts. This owes to the fact that Haar detail 
coefficients are insensitive to polynomial trends in the signal 
up to first order. Subtraction of a flat or linear background 
from a light curve is an operation under which the wavelet 
transform is invariant (as are Fourier transforms) apart from 
the mean signal coefficient. 

The GRB light curves show power at various variablity 
time scales. Most often, there is a region of the logscale di- 
agram (log-power verses log-varibility time) with a single 
slope, indicating scaling in the power over those variability- 
times, and a flat region at the shortest variability times, indi- 
cating the presence of white-noise. Some of this white-noise 
may be intrinsic to the GRB. Some may be attributed to in- 
strumental noise and to background emissions from sources 
not including the GRB in question. We therefore express the 



variability of the burst, Pj UTBt , at time scales j as compris- 
ing of individual variances: a scaling component, / g scahn s- a n 
intrinsic noise component, p nolsc - an d a background com- 
ponent, ^background ^ The variability f the burst can then 
be described as a linear combination of the component vari- 
ances so long as the components have vanishing covariances. 
In this event we write, 







^scaling ^ ^noise ^ ^background 
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The minimum variability time scale, rg, is identified from a 
logscale diagram by the octave, ^'intersection , of the intersec- 
tion of the flat intrinsic noise domain, /3" olsc , with the sloped 
scaling domain, , 



Tg = T bin x 2 



.7 intersection 



(8) 



In practice, the octave at which the intersection occurs is de- 
termined by equating the polynomial fits to the flat intrinsic 
noise domain and the sloped scaling domain and solving for 
^intersection- The uncertainty in t$ is determined by propa- 
gating the uncertainty in the parameters from the fits to 
the j3j which in turn follow from a bootstrap procedure de- 
scribed in Sec. 14.1.21 and Sec. 14.1.31 It is at this time scale, 
t q , that a structured physical process appears to give way 
to one that is stochastic and unstructured. Clearly one seeks 
to remove ^background £ rom [j] £ Q arr j ve a t the cleanest 

possible signal, 



^background ^scaling ^ ^noise 



(0) 



^burst ^iclean /iburst 

Pj —t Pj = Pj r-j r-j 

In order to estimate the variance of the background dur- 
ing the burst, we will assume that the variance obtained 
from a preburst portion of the light curve can serve as an 
acceptable surrogate for the background variance. That is, 
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preburst 
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background 



(10) 



and then the background is removed from the signal accord- 
ing to the relation, 



0? 



(11) 



A simple algebraic manipulation of Eq. [TT] gives a form, 



lo g2 (^ Iean )=log 2 (^ urBt ) + log 



oprcburst \ 

LA I (12) 

^burst j V ' 



For long GRBs, the preburst is defined relative to a s 
trigger time as T-20 s to T-5 s and for short GRBs the 
preburst is defined to be from T-15 s to T-l s. Here T is the 
trigger time of the burst. 

4-1.1 Statistical Uncertainties 

We have considered the statistical uncertainties in the light 
curve by a typical bootstrap approach in which the square 
root of the number of counts per bin is used to generate 
an additive poisson noise. A new poisson noise is considered 
for each iteration through the bootstrap process. More sig- 
nificant contributions to the uncertaintaies are discussed in 
SecgXHand Sec. l4~j~31 

4-1.2 Circular Permutation 

Spurious artifacts due to inc idental symmetries resulting 
from accidental misalignment (|Percivalll2o6"o1 ; ICoifmadll995f ) 
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of light curves with wavelet basis functions are minimized 
by circularly shifting the light curve against the basis func- 
tions. C ircular shifting is a form of translation invariant de- 
noismg (|Coifmadll995h . It is possible a shift will introduce 
additional artifacts by moving a different symmetry into a 
susceptible location. Thus, our approach is to circulate the 
signal through all possible values, or at least a representative 
sampling, and then take an average over the cases which do 
not show spurious correlations. 



4-1.3 Reverse-Tail Concatenation 

Both discrete Fourier and discrete wavelet transformations 
imply an overall periodicity equal to the full time-range of 
the input data. This can be interpreted to mean that for a 
series of N elements, {Xq,Xi . ..Xjv-i} then Xq is made 
a surrogate for Xn and X\ is made a surrogate for Xjf+i, 
and so forth. This assumption may lead to trouble if Xq 
is much different from Xm-i- In this case, artificially large 
variances may be computed. Reverse-tail concatenation min- 
imizes this problem by making a copy of the series which is 
then reversed and concatenated onto the end of the original 
series resulting in a new series with a length twice that of 
the original. Instead of matching boundary conditions like, 



Pseudo Light Curve 



Pseudo Logscale Diagram 



Xq, Xi, . . . , Xjv-i, A'o, 
we match boundaries as, 
Xo, Xi, . . . Xn-i, Xn-i, 



, X\ , Xq . 



(13) 



(14) 



Note that the series length has thus artificially been in- 
creased to 2N by reversing and doubling of the original se- 
ries. Consequently, the wavelet variances at the largest scale 
in a logscale diagram reflect this redundancy. This is the rea- 
son the wavelet variances at the largest scale are excluded 
from least-squares fits of the scaling region. 

Another difficulty in wavelet expansions is that the 
initialization procedure of the multi-resolution algorithm 
may pollute the detail coefficients at the finest scale (see 
IStrang and Nguvenlll997l; lAbrv et alJ|2003l ). For this reason 
we follow the advice of lAbrv et al. ( 20031 ) and discard the 
detail coefficients at the finest scale. 



4.2 Simulation 

The efficacy of this background subtraction method was 
tested using simulated data in the form of fractional Brown- 
ian motion ( fBm) time series, which were first discussed by 
iMandelbrotl | 19681 ). The fBms were combined with Poisson 
noise. The numerical computing environment MATLAB was 
used to produce 1000 realizations of fBms with scaling pa- 
rameter a randomly chosen from the range 1.0 ^ a ^ 2.0 by 
using a uniform random number generator. The fBms were 
combined with a Poisson noise with variance, Ai, randomly 
chosen from the range 0.25 ^ Ai ^ 1.0 by using a uniform 
random number generator. This Poisson noise is regarded as 
intrinsic to the GRB. The fBms and the Poisson noises thus 
produced were combined and stored as ideal light curves, 
shown in green on the left panel of Fig. [2] The logscale 
diagrams obtained from the ideal light curves are used to 
compare against the results of the subtraction procedure. 
Another set of Poisson noise with variances, A2, were 
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Figure 2. In the left hand panel are simulated light curves and 
noise processes: an ideal fBm process (green) and a background 
Poisson noise (black). The sum of the fBm and Poisson processes 
is shown in red and is labeled observed. The observed light curve 
(red) is the sum of the fBm and Poisson noise. The right hand 
panel shows the results of the background subtraction procedure. 
Red and black points show the logscale diagrams corresponding to 
the observed light curve and background, respectively. The green 
data shows the logscale diagram for the ideal light curve and 
the blue data are the logscale with background removed. The 
agreement between green and blue data demonstrates the merit 
of the background removal procedure. 



randomly chosen from the range 1.0 ^ A2 ^5 4.0 by using a 
uniform random number generator. These noise signals are 
interpreted as external background meant to be removed by 
the subtraction procedure and are shown in black on the left 
panel of Fig. [5] 

The ideal light curves were then combined with exter- 
nal backgrounds resulting in pseudo- observed light curves 
shown in red on the left panel of Fig. [2] The pseudo-observed 
light curves and the external background noises were trans- 
formed into wavelet coefficients and wavelet variances were 
computed according to Eq. [5] The variances of the pseudo- 
observed light curve (labeled actual) and the background 
are plotted in the right panel of Fig. [2] in red and black, re- 
spectively. The background was subtracted from the pseudo- 
observed light curve as detailed in Eq. \T%\ and the resulting 
corrected variances are plotted in blue in the right panel 
of Fig. [2] The corrected variances are to be compared to 
the variances of the ideal light curve which are plotted in 
green. Notice the close agreement between the background- 
subtracted variances and the variances of the ideal light 
curve. The pseudo data shown in Fig. [5] is one realization 
out of the 1000 that were generated but is representative of 
the set. 

In summary, 1000 simulated light curves were gener- 
ated and background noise was added. The light curves 
with background noise were then denoised using the same 
algorithm applied to actual GRB data in which preburst 
data were used as a surrogate for background. The simu- 
lated background subtracted variances were then compared 
to the variances of the ideal light curves, i.e., light curves 
without external background noise. The results indicate that 
the background subtraction method is robust and gives con- 
fidence that external background noise can be subtracted 
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Figure 3. A histogram of minimum variability time scales, in 
the observer frame, for long and short GRBs. It is clear that the 
distribution of long GRBs is displaced from the distribution of 
short GRBs. 



from the GRB light curves with the assumption that pre- 
burst data can serve as a surrogate for background noise. 



4.3 Selection Criteria 

We analyzed 122 GRBs (61 long and 61 short) listed in the 
Fermi GBM-Burst Catalog (jPaciesas et al.ll2012t ) for the first 
two years of the GBM mission. As discussed in Sec. 14.21 the 
signal-to-background ratio is a factor to be considered in 
recovering the intrinsic light curve (see Eq. ll2[) . We required 
the following condition on the ratio of variances, 

^prcburst 

^ < 0.75, 



ft 



(15) 



at all octaves, j. In addition, we also required that the first 
order polynomial fits to the noise region and to the scaling 
region each had a x 2 /d.f. that was less than 2. This reduced 
the sample to 14 short GRBs (Tab. [lj and 46 long GRBs 
(Tab. [2| for a total of 60 and it is these GRBs which are 
used to create Figs. and(7] For boosting into the source 
frame (Figs. [5] and [6]) a known z is obviously required and 
this cut further reduced the data set to 2 short GRBs and 16 
long GRBs for a total of 18 GRBs considered in the source 
frame (see Tab. [3]). 



5 RESULTS AND DISCUSSION 

For a large sample of short and long GBM bursts, we have 
used a technique based on wavelets to determine the min- 
imum time scale (rp) at which scaling processes dominate 
over random noise processes. The rp is the intersection of the 
scaling region (red-noise) of the spectrum in the logscale di- 
agram with that of the fiat portion representing the (white- 
noise) random noise component. This transition time scale is 
the shortest resolvable variability time for physical processes 
intrinsic to the GRB. Histograms of the extracted Tp values 
for long and short GRBs are shown in Fig. [3] We make two 
observations regarding these histograms: (1) There is a clear 
temporal offset in the extracted mean rp values for long and 
short GRBs. We believe this is the first clear demonstration 



Figure 4. Minimum variability time scale versus Tgo in the Ob- 
server frame. 
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Figure 5. Minimum variability time scale versus Tgn in Source 
frame. The correction for time dilation shortens Tgo and decreases 
the minimum variability time scale of each burst. 
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Figure 6. Minimum variability time scale versus Tgn with symbol 
size determined by luminosity (larger symbols for higher luminos- 
ity). No obvious relation between minimum variability time scale 
and luminosity is apparent. Sec Fig.[5]for error bars. 



6 MacLachlan et al. 



T90 v s T 90 /Tp (Observer Frame) 



10000 



1000 



100 



Long GRBs 
Short GRBs 




10 



0.1 



10 

T90 [sec] 



100 



1000 



Figure 7. The ratio of duration-to-minimum variability time 
scale (T90 /tr ) versus T90 . 



of this temporal difference. IWalker et al.l (|200Cj ). who stud- 
ied the temporal variability of long and short bursts using 
the BATSE data set did not report a systematic difference 
between the two types of bursts. (2) The two histograms are 
quite broad and very similar in dispersion. While the differ- 
ence in the mean rp is understandable (a point we discuss 
further elsewhere) the similarity of the dispersion is some- 
what surprising since the progenitors and the environment 
for the two types of bursts are presumably very different. 
The comparison is qualitative at best however because the 
rp scale has not been corrected for redshift (2), an effect that 
impacts the long bursts more than the short bursts. In pass- 
ing we note that the dispersion of the rp hist ogram (for long 
bursts ) is in agreement with the results of lUkwatta et all 
|201ll ) who performed a power density spectral analysis of 
a large sample of Swift long GRBs. In that work the au- 
thors extracted threshold frequencies and related them to a 
variability scale. 

In Fig. fj] we show a log-log plot of rp versus Tgn (the 
duration of the bursts); long GRBs are indicated by circles, 
the short ones by squares and both time scales are with re- 
spect to the observer frame. As in the histograms above, 
the fact that short GRBs, in general, tend to have smaller 
rp values compared to long GRBs, is evident in this figure. 
Also shown in the figure (as a dash line) is the trajectory 
of rp equal to T90. As we expect, no long GRBs exhibit a 
rp longer than T90 although interestingly a few short GRBs 
of extremely short duration appear to be approaching the 
limit of equality. In addition to establishing a characteristic 
time scale for short and long bursts, this figure also hints at 
a positive correlation between this time and the duration of 
bursts. We note that the rp scale spans approximately two 
decades for both sets of GRBs and that the two groups are 
fairly well clustered in the Tp-Tgo plane. A closer examina- 
tion of the two groups, however, indicates that a correlation 
between rp and T90, if present, is marginal at best. This is 
certainly true for the short-GRB group, especially given the 
large uncertainties in the Tgns for these bursts. The situation 
for the long-burst group on the other hand is not immedi- 
ately clear. In order to explore this further we cast the rp 
and the Tgn time scales into the source frame by applying 
the appropriate (1 + z) factor to the GRBs for which the z 



is known. Unfortunately the z is not available for the ma- 
jority of the short GRBs but we note that the correction 
is the same for both axes and is, to first order, small for 
the short GRBs since the mean z for this group is < 0.8. 
The corrected results for long-GRBs are shown as a log-log 
plot in Fig. [5] We see from this figure (and Fig. [4} the ap- 
pearance of a very intriguing feature: A plateau region in 
which the rp is essentially independent of T90 and a scaling 
region in which the rp appears to increase with Tgn, with 
the transition occurring around T90 on the order of a few 
seconds. 

If one assumes a positive correlation between luminosity 
and variability as suggested by a number of authors, then 
one might expect smaller rp values for higher luminosity 
bursts compared to those of lower luminosity. To investigate 
this, the data (in Fig. [5} are re-plotted in Fig. [6] in which 
the size of each datum symbol has been modulated by the 
gamma-ray luminosity of the burst, i.e., a large symbol im- 
plies a high luminosity and a small symbol a low luminosity. 
We see from Fig. [5] that no obvious connection between rp 
and luminosity is evident. 

Under the assumption that the rp is a measure pro- 
portional to the smallest causally-connected structure as- 
sociated with a GRB light curve, it is then possible to 
interpret the scaling trend in terms of the internal shock 
model in which the basic units of emission are assumed to 
be pulses that are produced via the collision of relativistic 
shells emitted by the central engine. Indeed, we note that 
iQuilligan et al. I (20021 ') in their study of the brightest BATSE 
bursts with Tgn > 2 sec explicitly identified and fitted dis- 
tinct pulses and demonstrated a strong positive correlation 
between the number of pulses and the duration of the burst . 
More recent studies dBhat fc Guiriedl201ll;lBhat et alj|2012l ; 
lHakkila fc Cumbed 120081 ; lHakkila fc Preecd l201lf ) provide 
further evidence for the pulse paradigm view of the prompt 
emission in GRBs. In our work we have not relied on iden- 
tifying distinct pulses but instead have used the multi- 
resolution capacity of the wavelet technique to resolve the 
smallest temporal scale present in the prompt emission. If 
the smallest temporal scale is made from pulse emissions 
from the smallest structures, then we can get a measure 
of the number of pulses in a given burst through the ratio 
Tgo/rp . In the simple model in which a pulse is produced 
every time two shells collide then the ratio, Tgo/rp, should 
show a correlation with the duration of the burst. A plot of 
this ratio versus Tgn is shown for a sample of short and long 
bursts in Fig. [7] The correlation is apparent. 

It is now widely accepted that the progenitors for the 
two classes of GRBs are quite distinct i.e., the merger of 
compact objects in the case of short GRBs and the collapse 
of rapidly rotating massive stars in the case of long GRBs. 
Formation of an accretion disk in the two cases is posed 
in a number of models but important factors such as the 
size of the disk, the mass of the disk, the strength of the 
magnetic field, in addition to the magnitude of the accre- 
tion rate during the prompt phase, remain largely uncertain. 
With contributions from intrinsic variability of the central 
engine or nearby shock-wave interactions within a jet, we 
should not be surprised to observe a systematic difference 
in the extracted variability time scales for long and short 
bursts, since the progenitors have different spatial scales. 
Knowing the variability timescales, we can estimate the size 
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of an assumed emission region. From Fig. [4] we note that 
the smallest temporal-variability scale for the short bursts 
is approximately 3 ms and that for the long bursts is ap- 
proximately 30 ms: These times translate to emission scales 
of approximately 10 s and 10 9 cm respectively. Our variabil- 
ity times and size scales are g enerally consistent with the 
findings of lWalker et al.l (|2000l ) although these authors also 
reported observing time scales as small as few microseconds. 
We find no evidence for variability times as low as a few mi- 
croseconds. 

iMorsonv etafl (|2010h modeled the behavior of a jet 
propagating through the progenitor and the surrounding 
circumstellar material and showed that the resulting light 
curves exhibited both short-term and longterm variability. 
They attribute the long-term variability, at the scale of few 
seconds, to the interaction of the jet with the progenitor. 
The short-term scale, at the level of milliseconds, they at- 
tribute to the variat ion in the activity of the central engine 
itself. Alternatively, IZhang fc Yanl |201ll ) consider a model 
in which the prompt emission is the result of a magnetically 
powered outflow which is self-interacting and triggers rapid 
turbulent reconnections to power the observed GRBs. This 
model also predicts two variability co mponents but interest - 
ingly and in contrast to the findings of lMorsonv et all (|201Cj ) 
, it is the slow component that is associated with the activ- 
ity of the central engine, and the fast component is linked to 
relativistic magnetic turbulence. While we are not in a posi- 
tion to distinguish between these two models it is intriguing 
nonetheless to note (see Fig. |4| that indeed there do appear 
to be two distinct time domains for the rp: a plateau region 
dominated primarily by short bursts although it includes 
some long bursts too, and a scaling region (i.e., a hint of a 
correlation between rp and T90) that is comprised solely of 
long bursts. In addition, we observe that the time scale in 
the plateau region is the order of milliseconds whereas that 
for the scaling region is approaching seconds. 

There is considerable dispersion in the extracted Tp. 
The variation is evident for both short and long-duration 
GRBs. The main cause of this dispersion is not fully un- 
derstood but one factor that may play a sig nificant role is 
angular momentum. As lLindner et al. |2010l ) note, the ba- 
sic features of the prompt emission can be understood in 
terms of accretion that results via a simple ballistic infall of 
material from a rapidly rotating progenitor. Material with 
low angular momentum will radially accrete across the event 
horizon whereas the material with sufficient angular momen- 
tum will tend to circularize outside the innermost stable cir- 
cular orbit and form an accretion disk. Simulations that go 
beyo nd the simple radial infall model (|Lindner et al.l [2010l . 
l201lf ) suggest that the formation of the disk leads to an ac- 
cretion shock that traverses outwards through the infalling 
material. If one assumes that the initiation of such an ac- 
cretion shock and the subsequent emission of the prompt 
gamma-rays are associated with a particular time scale, the 
variability of this scale then (as given by the dispersion in 
rp for example) may reflect the different dynamics (initial 
angular momentum and the mass of the black hole) of each 
GRB in our sample. In the case of long GRBs, the mass of 
the central black hole can vary by an order of magnitude thus 
potentially explaining a large part of the dispersion seen in 
the rp. However a similar dispersion for short bursts is dif- 
ficult to reconcile using the same arguments since the mass 



range for the central black hole in standard merger models 
(at least for NS-NS mergers) is expected to be significantly 
smaller. 



6 CONCLUSIONS 

We have studied the temporal properties of a sample of 
prompt-emission light curves for short and long-duration 
GRBs detected by the Fermi/GBM mission. By using a tech- 
nique based on wavelets we have extracted the variability 
timescales for these bursts. Our main results are summa- 
rized as follows: 

a) Both short and long-duration bursts indicate a tem- 
poral variability at the level of a few milliseconds. Variability 
of this order appears to be a common feature of GRBs. Thi s 
finding is consistent with the work of I Walker et al.l (|200(t ). 
However, unlike these authors we do not find evidence of 
variability at a time scale of few microseconds. 

b) In general the short-duration bursts have a variabil- 
ity time scale that is significantly shorter than long-duration 
bursts. In addition, the Tp values seem not to depend in 
any obvious way on the luminosity of the bursts. The dis- 
persion over different GRBs in the extracted time scale for 
short-duration bursts is an order of magnitude within the 
smallest variability time, that time being approximately 3 
milliseconds. The dispersion for the long-duration bursts is 
somewhat larger. The origin of the dispersion in either case 
is not known, although we should expect that the size of the 
initial angular momentum and the mass of the system play 
significant roles. 

c) The ratio of Tgo/rp appears to be positively corre- 
lated with the minimum variability time scale. This suggests 
further support for the pulse paradigm view of the prompt 
emission as being the result of shell collisions. In this re- 
spect, the minimum variability time scale is likely related to 
key pulse parameters such as risetimes and/or widths. 

d) For short-duration bursts, the variability parameter 
773 shows negligible dependence on the duration of the bursts 
(characterized by T90). In contrast, the long-duration bursts 
indicate evidence for two variability time scales: a plateau 
region (at the shortest time scale) which is essentially inde- 
pendent of burst duration and a scaling region (at the higher 
time scale) that shows a positive correlation with burst dura- 
tion. The transition between the two regions occurs around 
T90 on the order of a few seconds in the source frame. 
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Table 1. Short GRBs (Observer Frame). 



GRB 


Tgo [sec] 


8T< )a [sec] 


n 


: [sec] 


6r t 


I [sec] 


Sr f 


t i sec \ 


UoU / zoy lo 





.192 





,345 


U 


UoU / 


n 

u. 




U 


UOlU 


081012045 


1 


216 


1 


,748 





,0052 





0024 
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1 
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231 





0258 
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,0165 


081105614 


1 


,280 


1 


,368 





0306 





0147 





0282 


081107321 


1 


,664 





,234 





,0504 





0129 





,0173 


081216531 





,768 





,429 





,0138 





0037 
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090108020 





,704 





,143 





,0241 





0064 





0088 


090206620 
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,143 





0143 





0063 
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,0005 
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1 


,664 
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0040 
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090429753 





,640 





,466 





0285 





0115 





0193 


090510016 





,960 





138 





,0049 





0009 





0011 


100117879 





,256 





,834 





0331 





0122 





,0192 
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Table 2. Long GRBs (Observer Frame). 



GRB 


T 90 [sec] 
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1 [ sec ] 
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020790. 
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c;q QftQ 
Do.oOy 
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2 
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CI 

u. 


D9Q9 
UZoZ 
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1 H7 /I K*7 
1U / .40 f 
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1 
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CI 

u. 
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U 
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1 
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u. 
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0. 
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u 
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CI 

U. 
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CI 
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U 
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o 
u. 
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0. 
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U 
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O 
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Ci 

u. 
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1 / 4o 


CI 

U. 
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o 
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/1 1 Q/l c; 
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CI 

u. 
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ci 
u 
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Uol IUIOOZ 


o.ZOO 


U 


ooy 


n 
u 


uy4o 
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U 


OQOO 

UoUZ 


o 
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uoiizo iyu 


O 920 

y.zou 
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n 
u 
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ZloZ 


(J. 


0504 


ci 
u 


Oft f^ft 
UDOO 
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(^9 fi^7 


7 
( 


qi o 


C| 

u. 
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uy iz 


CI 

U. 


oooo 
uzyz 


o 
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O.OOo 
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n 
U 


Uoiy 


0, 


0043 


ci 
u 
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UUOU 
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U 
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ci 

u. 
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CI 

U. 
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Uo41 


o 
u 


Uo41 


021 99990/1 
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2 
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o 
u. 
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0. 
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ci 
u 


07Q9 
U i OZ 
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10.44O 


1 


.159 


o 
u. 
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0. 
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ci 
u 
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U4ol 
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uyuiuzizz 
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ZD.0Z4 
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CI 

U. 
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Uo4( 


0. 
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ci 
u 
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U1D4 
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07Q 
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1 
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CI 

U. 
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U / OO 


0. 
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u 


0990 
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uyuzuzo4 ( 
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u. 
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0. 
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u 
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loo. 1 / U 


1 
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u 


ioyo 
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0436 


n 
u 


uoyy 
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1 
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CI 

u. 
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0. 
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u 
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14.00U 


1 


.086 


CI 

u. 
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0. 
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u 
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.264 


o 
u. 
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UZ4y 


0. 
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U 


OOQft 
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2 


.450 


o 
u. 
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0. 
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U 
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1 


.028 


o 
u 
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u 
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2 
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u 
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i 

i 
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1 
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u. 
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0. 
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U 
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U4ZO 
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1 IZ.OoO 


1 


,086 


n 

u. 
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ZDOl 





0536 


n 
U 


Oft7Q 
UU I O 


090620400 


13.568 





.724 


0. 


1667 


0. 


0422 





.0564 


090626189 


48.897 


2 


.828 


0. 


0498 


0. 


0078 





.0093 


090718762 


23.744 





.802 


0. 


1621 


0. 


0482 





0686 


090809978 


11.008 





,320 


0. 


2436 


0. 


0515 





0652 


090810659 


123.458 


1 


.747 





7319 





3027 





.5161 


090829672 


67.585 


2 


.896 


0. 


0678 


0. 


0141 





.0177 


090831317 


39.424 





,572 


0. 


0266 


0. 


0103 





.0169 


090902462 


19.328 





,286 


0. 


0223 


0. 


0026 





0029 


090926181 


13.760 





.286 


0. 


0435 


0. 


0061 





.0070 


091003191 


20.224 





.362 


0. 


0300 


0. 


0051 





.0062 


091127976 


8.701 





,571 





0395 





0059 





0069 


091208410 


12.480 


5 


,018 


0. 


0621 


0. 


0180 





.0254 


100414097 


26.497 


2 


.073 


0. 


0418 


0. 


0074 





.0090 
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Table 3. Long and Short GRBs (Tgo and rp in Observer Frame). Luminosities are taken from references given in footnotes. 



GRB 




z 


T go [sec] 


<5T 90 [sec] 




[sec] 


Sr. 




5rj" [sec] 


L lso [ergs/s] 


SL iso [ er S S / S l 


lo g SLf o [ergs/s] 


0808019/ 2 


2 


204 


24.704 


1.460 





4306 





1336 


0.1937 


3.5810 


FT O O "1 o n 1 

5.8210 


7.85-10 


08081U549 


3 


350 


107.457 


15.413 





1353 





0648 


0.1243 


9.5911) 


i o o i o n 2 


1 o o 1 n52 


f\c*](\t\ 1 f c\r\c\ 

080916009 


1 


350 


b2.9l" 





810 





2266 





0630 


r\ r\o to 

0.0872 


3 1 o A 1 Ao4 

°1. 04-10 


8.79-10 


O f ~7f\ "I 0^2 

8. 79 10 


AOI oooon/i 


2 


770 


18.880 


2 


318 





1956 





0533 


0.0732 


1 1 oc i r\53 


7-10 


6-10 51 


090102122 


1 


547 


26.624 





810 





0347 





0111 


0.0164 


^71 1 O Pi 2 

8.71-10 


5.6-10 51 


5.6-10 51 


090323002 


3 


570 


135.170 


1 


448 





1598 





0436 


0.0599 


1 6.87-10 53 


6.55-10 53 


4.45-10 52 


090328401 





736 


61.697 


1 


810 





0682 





0139 


0.0175 


4 1.79-10 52 


1.42-10 51 


1.11-10 51 


090424592 





544 


14.144 





264 





0249 





0031 


0.0036 


5 1.62-10 52 


4-10 50 


5-10 50 


090510016 





903 


0.960 





138 





0049 





0009 


0.0011 


3 1.78-10 53 


1.2-10 51 


1.2-10 51 


090516353 


4 


100 


123.074 


2 


896 





7992 





5686 


1.9711 


6 8.17-10 52 


2.85-10 52 


6.1-10 51 


090618353 





540 


112.386 


1 


086 





2631 





0536 


0.0673 


5 8.47-10 51 


1.1710 51 


3.4-10 50 


090902462 


1 


822 


19.328 





286 





0223 





0026 


0.0029 


3 5.89-10 53 


9.71-10 51 


9.71-10 51 


090926181 


2 


106 


13.760 





286 





0435 





0061 


0.0070 


3 7.40-10 53 


1.45-10 52 


1.45-10 52 


091003191 





897 


20.224 





362 





0300 





0051 


0.0062 


^^■lO 52 


3.71-10 51 


6.55 10 51 


091127976 





490 


8.701 





571 





0395 





0059 


0.0069 


7 3.70-10 51 


1.38-10 50 


1.06-10 50 


091208410 


1 


063 


12.480 


5 


018 





0621 





0180 


0.0254 


^^■lO 52 


1.48-10 51 


3.45-10 51 


100117879 





920 


0.256 





834 





0331 





0122 


0.0192 


^.es-io 52 


5.0110 51 


1.08-10 52 


100414097 


1 


368 


26.497 


2 


073 





0418 





0074 


0.0090 


8 1.00-10 53 


1.58-10 52 


7.6-10 51 



1 Nava et al. 2011 (A&A 530, A21 (2011)) 

2 GCN #8100 

3 Ghirlanda et al. 2011 l larXiv:1107.4096l l 

4 GCN #9057 

5 Ukwatta, T. N., et al. 2010, ApJ., 711, 1073 

6 GCN #9415 

7 GCN #10204 

8 GCN #10595 
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